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In this work we review some features of topological defects in field theory models for real scalar 
fields. We investigate topological defects in models involving one and two or more real scalar fields. 
In models involving a single field we examine two different subclasses of models, which support 
one or more topological defects. In models involving two or more real scalar fields, we explore the 
presence of defects that live inside topological defects, and junctions and networks of defects. In 
the case of junctions of defects we investigte structures that simulate nanotubes and fulerenes. Our 
investigations may also be used to describe nonlinear properties of polymers, Langmuir films and 
optical solitons in fibers. 



I. INTRODUCTION 

Research on topological defects in Field Theory was 
iniciated almost three decades ago, and some of the main 
investigations may be found for instance in Refs. 
In the case of topological defects that appears in models 
described by real scalar fields in (1,1) space-time dimen- 
sions, they are usually named kinks, and are classical 
static solutions of the equations of motion. Their topo- 
logical behavior is related to the asymptotic form of the 
field configurations, which has to differ in both the pos- 
itive and negative space directions. To ensure that the 
classical solutions have finite energy, one requires that the 
asymptotic behavior of the solutions is identified with 
minima of the potential that defines the system under 
consideration, so in general the potential has to include 
at least two distinct minima for the system to support 
topological solutions. 

We can investigate real scalar fields in (3, 1) space- 
time dimensions, and now the topological solutions are 
named domain walls. These domain walls are bidimen- 
sional structures that carry surface tension, which is iden- 
tified with the energy of the classical solutions that spring 
in (1,1) space-time dimensions. The domain wall struc- 
tures are supposed to play a role in applications to several 
different contexts, ranging from the low energy scale of 
condensed matter ^| || up to the high energy scale re- 
quired in the physics of elementary particles, fields and 
cosmology 

There are at least three classes of models that support 
kinks or domain walls. In the first class one deals with 
a single real scalar field, and the topological solutions 
are structureless. Examples of this are the sine-Gordon 
and (p 4 models. In the second class of models one also 
deals with a single real scalar field, but now the systems 
comprise at least two distinct domain walls. An example 
of this is the double sine-Gordon model, which has been 
investigated for instance in Refs. |H-ftj]. In the third 
class of models we deal with systems defined by two real 
scalar fields, and now one opens two new possibilities: 
domain walls that admit internal structure p^pp|, and 
junctions of domain walls, which appear in models of two 



fields when the potential contains non-colinear minima, 
as recently investigated for instance in Refs. ]2~H-|35|]. 

There are other motivations to investigate domain 
walls in models of field theory, one of them being related 
to the fact that the low energy world volume dynamics of 
branes in string and M theor y m ay be described by stan- 
dard models in field theory [p6|-p8[. Besides, one knows 
that field theory models of scalar fields may also be used 
to investigate properties of quasi-linear polymeric chains, 
as for instance in the applications of Refs. [|39|-fl2||, to de- 
scribe solitary waves in ferroelectric crystals, the pres- 
ence of twistons in polyethylene, and solitons in Lang- 
muir films. 

In the present work, in Sec. || we review some known 
facts about models described by real scalar fields. The 
investigations follow in Sees. fTD and IV, where we search 



for topological structures that generate kinks and walls, 
in models involving a single real scalar field, and two or 
more fields, respectively. 



II. GENERAL CONSIDERATIONS 

In this work we are interested in field theory models 
that describe real scalar fields and support topological 
solutions of the Bogomol'nyi-Prasad-Sommerfield (BPS) 
type |^3|,0 • In the case of a single real scalar field <p, we 
consider the Lagrange density 
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Here V{<p) is the potential, which identifies the particu- 
lar model under consideration. We write the potential 
in the form V{cf>) = (1/2) W|, where W = W{<p) is a 
smooth function of the field (j>, and W$ = dW/d<p. In a 
supersymmetric theory W is the superpotential, and this 
is the way we name W in this work. 

The equation of motion for <p = <p{x, t) has the general 
form 
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and for static solutions we get 

d£ = W + W ** 



(3) 



It was recently shown in Ref. |15| that this equation of 
motion is equivalent to the first order equations 



(4) 



if one is searching for solutions that obey the boundary 
conditions lim^^-oo cf>(x) = 4>i and ]im x —,- 00 (d(f)/dx) = 
0, where 4>i is one among the several vacua {0i,</>27 ■■■} 
of the system. In this case the topological solutions are 
BPS (+) and anti-BPS (-) solutions. Their energies get 
minimized to the value ft = | AWy |, where AWjj = Wj — 
Wj, with Wi standing for Wfa). The BPS and anti-BPS 
solutions are defined by two vacuum states belonging to 
the set of minima that identify the several topological 
sectors of the model. 

In the case of two real scalar fields and x the potential 
is written in terms of the superpotential, in a way such 
that V((f>,x) = (1/2) Wl + (1/2) W%. The equations of 
motion for static fields are 
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which are solved by the first order equations 



f=±W, £ 
dx dx 
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Solutions to these first order equations are BPS (+) and 
anti-BPS (-) states. They solve the equations of mo- 
tion, and have energy minimized to ft = |AW y '| as 
in the case of a single field; here, however, AW 1 ^ = 
W((f>i, Xi) — W(4>j,Xj), since now we need a pair of num- 
bers (</>j, Xi) to represent each one of the vacuum states 
in the system of two fields. In the plane (4>,x) we m ay 
have minima that are non colinear, openning the pos- 
sibility for junctions of defects. In the case of two real 
scalar fields, we can find a family of first order equations 
that are equivalent to the pair of second order equations 
of motion, but this requires that the superpotential is 
harmonic, obeying + W xx = |45) . 



III. MODELS INVOLVING A SINGLE REAL 
SCALAR FIELD 

We now turn attention to kinks and domain walls. A 
well known example is given by the <j> model, defined by 
the potential V(<f>) — (1/2) ((f) 2 — l) 2 . Here we are using 
natural units, and dimensionless fields and coordinates. 



In this model the domain wall can be represented by the 
solution 4> s (x) = ±tanh(x). The above potential can be 
written with the superpotential W(4>) — <fi — 3 /3, and 
the domain wall is of the BPS or anti-BPS type. The 
wall tension corresponding to the BPS wall is t s — 4/3. 

We can also find structureless domain walls in other 
models, for instance in the 4> e model, which is described 
by the potential V(<j>) = (1 / '2) cf> 2 (4> 2 - I) 2 . Here we 
have W(4>) = (l/2)0 2 - (l/4)</> 4 , and the wall config- 
urations are also of the BPS type, and are given by 
4>l = (l/2)[l±tanh(ar)]. The wall tension is now t s = 1/4. 
This potential was investigated for instance in Ref. |4q ]. 

We can build another class of models where the do- 
main walls engender other features. The next class is yet 
described by a single field, but the systems may now sup- 
port two or more different walls. An interesting example 
of this is the double sine-Gordon model, which is defined 
by the potential 
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V r U) = \4r cos(6) + cos(20)l 

r + 1 



(8) 



where r is a parameter, real and positive. This potential 
is periodic, with period 27r; for simplicity in the following 
we consider the interval — 2tt < <f> < 2ir. The value r = 1 
distinguishes two regions, the region r£ (0,1) where the 
potential contains four minima, and the region r > 1, 
where the potential contains two minima. For r 6 (0, 1) 
the system supports two distinct wall configurations, the 
large wall and the small wall, which distinguish the two 
different barrier the model comprises in this case. The 
limits r — > and r — > oo lead us back to the sine-Gordon 
model. The double sine-Gordon model has been consid- 
ered in several distinct applications, as for instance in 
Ref. P [Tl|] , where one investigates magnetic solitons in 
superfluid 3 He, kink propagation in a model for poling in 
polyvinylidine fluoride, and properties related to the two 
different kinks that appear in such polymeric chain. 

To expose new features of the double sine-Gordon 
model we rewrite Eq. (||) in the form |47]] 
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where we have omitted an unimportant r-dependent con- 
stant. The model can be described by the superpotential 



W{4>) 



[sin(0) + r 4>] 



(10) 



For r in the interval r £ (0, 1) the minima of the potential 
are the singular points of the superpotential, dW/d<fi = 0. 
They are periodic, and for —2tt < <f> < 2ir there are four 
minima, at the points <j> = ±7r ± a(r), where a(r) — 
cos _1 (r). For r > 1 the minima are at <f> = ±7r, in the 
interval — 2ir < <j> < 2tt. A closer inspection shows that 
for < r < 1 the local maxima at ±7r and the minima 
±7r ± a(r) degenerate to the minima ±7r for r = 1, and 
remain there for r > 1. Thus, the parameter r induces 
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a transition in the behavior of the double sine-Gordon 
model. The value r = 1 is the critical value, since it is the 
point where the system changes behavior: for r E (0, 1) 
this model supports minima that desapear for r > 1. We 
illustrate the double sine-Gordon model in Fig. 1, where 
we depict the potential of Eq. (||) for r = 1/3,2/3 and 
for r = 1. 
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FIG. 1. The double sine-Gordon potential, depicted for 
r = 1/3,2/3 and 1 from above to below respectively, to il- 
lustrate how the behavior of the model changes with r. 

We get a better view of the phase transition in the dou- 
ble sine-Gordon model by examining the order parameter 
4>(r), which is given by ±7r ± a(r) for < r < 1, so it 
goes continuously to ±7r for r > 1. Also, the (squared) 



mass of the field can be obtained via the relation 



(11) 



where cf> is the corresponding minimum of the potential. 
For < r < 1 we get m 2 (r) — 4 — 4r, and for r > 1 we 
have m 2 (r) = 4(r — l)/(r + 1). We plot <j>{r) and m 2 (r) 
in Fig. 2. We see that m(r) vanishes in the limit r — > 1. 
These results indicate that r drives a second order phase 
transition, a transition where the system goes from the 
case of two distinct phases to another one, engendering 
a single phase. 




FIG. 2. Plots of 4>{r) (above) and m 2 (r) (below) for the 
double sine-Gordon potential, which illustrate how the be- 
havior of the model changes with r. 

We consider < r < 1. The energies of the BPS 
solutions are given as follows. For solutions connecting 
the minima — 7r + a(r) and n — a(r) the defect is large 
since it joins minima separated by a higher and wider 
barrier. We have 



i . i 7t — air) 

yl + r 
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In the case of the minima tt — a(r) and 7r + a(r) the defect 
is small and we get 



b dsG 
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a(r) 
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We notice that t 



l 

dsG 



limit r 



1 sends t 



tdsG + 47rr/Vl + T, and that the 
"dsG -» 2 ^2vr and i* sG -> 0, as ex- 
pected. For the BPS states we can write the solutions 
explicitly. For instance, for solutions that connect the 
minima — 7r + a(r) and tt — a(r) we get large kink solu- 
tions, which are of the form 



4>i(x) = ±2 tan 



h+r 

V t - 7 tanh (y^* x ) 
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For solutions that connect the minima n ± a(r) and the 
minima — 7r±a(r) we get small kink solutions. They are 
given by 



,(x) = ±7r — 2 tan 



\ tanh — r x) 

V 1 + r 



(15) 



The potential in Eq. (H) in the limit r — > goes to 



V (4>) = l + cos(20) 



(16) 



which leads us back to the sine-Gordon model. Thus, we 
can suppose r small and use V r (0) to explore the double 
sine-Gordon model as a model controlled by a small pa- 
rameter, in the vicinity of the sine-Gordon model. This 
feature may be of some use for investigations that fol- 
low the lines of Ref. |Q , and also in the case concerning 
the presence of internal modes of solitary waves, which 
seems to appear when one slightly modifies some inte- 
grable model - see for instance Ref. G9) . 



IV. MODELS INVOLVING TWO OR MORE 
REAL SCALAR FIELDS 

We now turn attention to the third class of models, 
which is described by two real scalar fields. In this case 
the domain walls may engender internal structure. This 
line of investigation follows as in Refs. 0-0 and we 
illustrate such possibility with the system defined by the 
potential 



V{4>,x) = \{<t> l -l? + \r 2 



r(l + 2r)0 2 x 2 
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where the parameter r ^ is real. This model was first 
investigated in Ref. ]50| , and can be used to modify the 
internal structure of domain walls - see Refs. |]l6| , pT| . 
This potential follows from the superpotential 



W((j ) ,x) = ^-^ 3 ~r<P X 2 



(18) 



and the system supports two-field solutions. An explicit 
solution is 



4>{x) — tanh(2ra;) 
x(x) — a(r) sech(2rx) 



(19) 
(20) 



with a 2 (r) = 1/r — 2. This is a BPS solution, and now 
the parameter r is restricted to the interval r e (0, 1/2). 
The limit r — > 1/2 lead us to the one-field solution 4>{x) = 
tanh(x) and x( x ) — 0- As we know, all the BPS solutions 
are linearly stable |p3| . 

The two-field solutions obey 



+ x 2 /a 2 = i 



(21) 



which describes an elliptic arc connecting the two min- 
ima (±1,0) of the corresponding potential in the (<p,x) 
plane. The one-field solutions represent standard domain 
walls, while the two-field solutions may be seen as do- 
main walls having internal structure: the vector (0, x) hi 
configuration space describes an straight line segment for 
the one-field solution, and an elliptic arc for the two-field 
solution, resembling light in the linearly and elliptically 
polarized cases, respectively. The same solutions appear 
in condensed matter, in the anysotropic XY model used 
to describe ferromagnetic transition in magnetic systems, 
and there they are named Ising and Bloch walls, respec- 
tively - see for instance Ref. ||, chapter 7. 

We recall that domain walls may be seen as seeds [||J|] 
for the formation of non-topological structures. This 
possibility appears in Refs. ]5^-pq], where the discrete 
symmetry is changed to an approximate symmetry, or in 
Ref. | p7| , with the discrete symmetry biased so that do- 
mains of distinct but degenerate vacua spring unequally. 
Domain walls that appear in models involving two real 
scalar fields present features that are not seen in the case 
of a single field. Thus, we are now investigating other 
extensions of the above model, defined by the superpo- 
tential of Eq. (|l8|). In the work in Ref. [p8| we are in- 
vestigating the possibility of changing the elliptic arc to 
a more general orbit, that may find applications in con- 
densed matter. In another work |39[ we are investigating 
other superpotentials, non-polynomial, in an effort to ex- 
tend the model studied in |60 6|] to the case of two or 
more fields. The model defined by the superpotential of 
Eq. ( [Tsl ) present other interesting features, recently ex- 
plored in Ref. f33| , which deserve further investigations. 

We now turn our attention to polynomial potentials 
that engenders the Z3 symmetry, and that supports sta- 
ble three-junctions that generate a regular hexagonal net- 
work of defects. Investigations on the presence of junc- 
tions |Q of defects in supersymmetric models have been 
given in Refs. pl|-p3t, and in Ref. |25| we have found the 
fourth-order polynomial potential 
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This potential does not allow for supersymmetric exten- 
sions. The equations of motion for static configurations 
that follows in this case are 
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The potential has three degenerate minima, at the points 
«i = (3/2) (1,0) andw 2 ,3 = (3/4) (-1, ±V3). These min- 
ima form an equilateral triangle, invariant under the Z% 
symmetry. The distance between the minima is (3/2)y/3. 

We can obtain the topological solutions explicitly. The 
easiest way to do this follows by first examining the sector 
that connects the vacua v-i and v^. This is so because in 
this case we set cj) — —3/4, searching for a strainght-line 
segment in the (<fi, x) plane. This is compatible with the 
Eq. (p3|), and reduces the other Eq. ( |24| ) to the form 
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This implies that the orbit connecting the vacua 1)2 and 
V3 is a straight line. It is such that, along the orbit the 
X held feels the potential A 2 [x 2 ~ (27/16)] 2 . This shows 
that the model reduces to a model of a single field, and 
the solution satisfies the first-order equation 
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The solution is 
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x(x) — — -Vstanh \ — A: 



(26) 
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The other solutions can be obtained by rotations obeying 
the Z3 symmetry of the model. 

The full set of solutions of the equations of motion are 
collected below. In the sector connecting the minima Vi 
and V3 they are 
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In the sector connecting the minima V\ and vi they are 
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In the sector connecting the minima V\ and V3 they are 

(32) 
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The label (±) is used to identify kink and antikink. All 
the solutions have the same energy, (9/4) y / 27/8|A|. 

We examine how the bosonic fields behave in the back- 
ground of the classical solutions. We do this by consid- 
ering fluctuations around the static solutions 4>{x) and 
x{%)- We use the equations of motion to see that the 
fluctuations depend on the potential 
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v x4> v xx 
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Evidently, after obtaining the derivatives we substitute 
the fields by their classical static values 4>(x) and xi x )- 
The model under consideration is defined by the potential 
l7|). In this case we use ( p8| ) and ( p9| ) to obtain two 
ecoupled equations for the fluctuations. The potentials 
of the corresponding Schrodinger-like equations are 



Un{x) 
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4-2 sech^ 
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(35) 
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The eigenvalues can be obtained explicitly: in the x di- 
rection we get Wq = and w* = (9/2)y / X 2 /2, and in 

the 4> direction we have wf, = (9/2) y/X 2 /2. This shows 
that the pair ( p8| ) and (^9|) is stable, and by symmetry 
we get that all the three topological solutions are stable 
solutions. 

The classical solutions present the nice property of hav- 
ing energy evenly distributed in their gradient (g) and 
potential (p) portions. In terms of energy density they 
are 



5(s) = P(*) = \[Y 



A 2 sech 4 




(37) 



To understand this feature we recall the calculation done 
explicitly in the sector with <j) — —3/4, constant. There 
the model is shown to reduce to a model of a single field, a 
model that supports BPS solutions. Within this context, 
the above solutions are very much like the non-BPS so- 
lutions that appear in supersymmetric systems p6| ]. We 
use this property and the topological current 



J a = e al3 d, 



(38) 



It obeys d a J a — 0, and it is also a vector in the (<fi,x) 
plane. For static configurations we have J l a J a — p t p, 
where p = p{cf>, x) is the charge density. This charge 
density allows writing p l p = e{x), where e(x) = g(x) + 
p(x) is the (total) energy density of the solution. We 
use this result and the notation ij, to identify the sector 
connecting the vacua (<fii,Xi) an d (^jiXj)i to show that 
for any two different sectors ij and jfc, i,j,k = 1,2,3 we 
get that 
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(Pij + pjkf {Ptj + Pjk) < phpij + PjkPjk 



(39) 



This condition shows that the three-junction is a process 
of fusion of defects that occurs exothcrmically, providing 
stability of junctions in the present model. This result 
is more general than the one in Ref . p4|| , which appears 
within the context of supersymmetry. Evidently, our re- 
sult also works for BPS and non-BPS solutions that ap- 
pears in supersymmetric models, with the property of 
having energy evenly distributed in their gradient and 
potential portions Q . 

We notice that the orbits corresponding to the stable 
defect solutions form an equilateral triangle in the ((f>, x) 
plane. This is so because the solutions are straight-line 
segments joining the three vacuum states in configuration 
space. They are degenerate in energy, and this allows 
associating to each defect the same tension 



(40) 



This makes iy < tjk +tk%,i,j,k = 1,2,3, and now the in- 
equality is strictly valid in this case, stabilizing the three- 
junction that appears in this model when one enlarges the 
space-time to three spatial dimensions. 

We consider the possibility of junctions in the plane, 
which may give rise to a planar network of defects. We 
work in (2, 1) space-time dimensions, in the plane (x,y). 
We identify the plane (x, y) with the space of configu- 
rations, the plane ((f), x)- We illustrate this situation by 
considering, for instance, the solutions we have already 
obtained. They are collected in Eqs. (||)-@ in (1,1) 
dimensions. In the planar case they change to 
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(43) 
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(45) 



|V3±^tanh(i^A(y-V3x)j (46) 



These planar defects are domain walls, and can be used 
to represent the three-junction in the limit of thin walls. 

The three-junction that appears in this ^-symmetric 
model allows building a network of defects, precisely in 



the form of a regular hexagonal network, as depicted in 
Fig. 3 in the thin wall approximation. In this network 
the tension associated to the defect is the tipical value 
of the energy in this tiling of the plane with a regular 
hexagonal network, which seems to be the most efficient 
way of tiling the plane. As we have shown, our model be- 
haves standardly in (3,1) dimensions. It supports stable 
three-junctions that generate a stable regular hexagonal 
network of defects. 




FIG. 3. A regular hexagonal network of defects, formed 
by three-junctions surrounded by domains representing the 
vacua vi — 1, V2 = 2, and v$ — 3. 

In Ref. |2^] the idea of nesting a network of defects 
inside a domain wall has been presented. A model that 
contains the basic mechanisms behind this idea was in- 
troduced in Ref. |H| . It is described by three real scalar 
fields a, (f>, and x, and is defined by the (dimensionless) 
potential, 



V(a,ct>,x) = 
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~ 4 

x 2 ) 2 



TO 1 — 



3X 2 



(47) 



Here r couples a to the pair of fields ((f), x). This poten- 
tial is polynomial, and contains up to the fourth order 
power in the fields. Thus, it behaves standardly in (3, 1) 
space-time dimensions. Also, it presents discrete Z2 x Z% 
symmetry. We set ((f), x) ~ * (0,0), to get the projection 
V(cr,0,0) -> V(a) = (2/3) (a 2 - 9/4) 2 . The projected 
potential presents Z 2 symmetry, and can be written with 
the superpotential W(a) = (2 v / 3/9)cr 3 - (3^3/2) a, in the 
form V = (l/2)(dW/da) 2 . The reduced model supports 
the explicit configurations <Jh(z) = ± (3/2) tanh(\/3z). 
The tension of the host wall is th = 3v3 = (3/2) nih, 
where rrih represents the mass of the elementary a me- 
son. Also, the width of the wall is such that lh ~ 1/V3- 
The potentials projected inside (a — > 0) and outside 
(a — > ±3/2) the host domain wall are Vi n ((f>,x) and 
V ou t(<f>, x)- Inside the wall we have 
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^n(0,x) = (0 2 +X 2 ) 2 -<M0 2 -3 X 2 ) 

This potential engenders the symmetry, and there are 
three global minima, at the points v" L — (3/2)(l,0) and 
v 2\ = (3/4)(— 1, ±\/3), which define an equilateral tri- 
angle. Outside the wall we get 

v out (<t>, X ) = (4> 2 + x 2 ? - <t>(<i> 2 - 3 x 2 ) 

+ \{r-l){<p 2 + X 2 ) (49) 

V ou t also engenders the symmetry, but now the min- 
ima depend on r. We can adjust r such that r > 9/8, 
which is the condition for the fields (f> and x to develop no 
non-zero vacuum expectation value outside the host do- 
main wall, ensuring that the model supports no domain 
defect outside the host domain wall. The restriction of 
considering quartic potentials forbids the possibility of 
describing the Z 3 portion of the model with the complex 
superpotential used in plfl . 

We investigate the masses of the elementary (f> and x 
mesons. Inside the wall they degenerate to the single 
value mm — 3\/3/2. Outside the wall, for r > 9/8 they 
also degenerate to a single value, m out (r) = 3a/ (r — l)/2, 
which depends on r. We see that m ou t{r = 4) = TOj n . 
Also, m out (r) > m in for r bigger than 4, and m out (r) < 
rriin for r in the interval (9/8,4). 

We study linear stability of the classical solutions 
(7 = <7h{z) and (4>, x) = (0,0). The fields (f> and \ van- 
ish classically, and their fluctuations (77™, £ n ) decouple. 
The procedure leads to two equations for the fluctuations, 
that degenerate to the single Schrodinger-like equation 

d 2 ib„(z) 9 ... . , , o , „ , 
I^ + 2 V{Z)MZ) = w ^ n ^ ( 5 °) 

Here V(z) = — 1 + r tanh 2 \/3z. This equation is of the 
modified Poschl- Teller type, and can be examined analyt- 
ically. The lowest eigenvalue is w% = (3/2) a/6 r + 1 — 6. 
There is instability for r < 5/2, showing that the host do- 
main wall with (cf), x) = (0, 0) is unstable and therefore 
relax to lower energy configurations, with (<fi,x) 7^ (0, 0) 
for r < 5/2. Inside the host domain wall the sigma 
field vanishes, and the model is governed by the potential 
Vi n ((j},x), which consequently may allow the presence of 
non-trivial (4>,x) configurations. The host domain wall 
entraps the system described by Vi n ((f>, x) f° r the param- 
eter r in the interval (9/8, 5/2). In this interval we have 
m-out < m-in, showing that it is not energetically favor- 
able for the elementary <f> and x mesons to live inside 
the wall for r e (9/8,5/2). The model automatically 
suppress backreactions of the 4> an d X mesons into the 
defects that may appear inside the host domain wall. 

In Ref. [£5j the potential inside the wall was shown to 
admit a network of domain walls, in the form of a hexago- 
nal array of domain walls. In the thin wall approximation 
the network may be represented by the solutions 



3 9 (1/27 \ 
01 = 8 + 8 tanh 1 2 V Y^ + ^J ( 51 ) 

3 3 ( 1 / 27 \ 

Xl = 8^ _ 8^ tanh ( 2VT (y + ^ x) ) (52) 

and by (4>kiXk)t obtained by rotating the pair (<f>i, xi) by 
2(fc— l)7r/3, for k = 2, 3. We identify the space (<f), x) with 
(x,y), so rotations in ((f), x) also rotates the plane (x, y) 
accordingly. The energy or tension of the individual de- 
fects in the network is given by, in the thin wall approx- 
imation, t n = (27/8)a/3/2 = (9/8) m in . In the nested 
network, the width of each defect obeys l n ~ a/8/27. 
This shows that lu/l n = 3/2 \/2, and so the host do- 
main wall is slightly thicker than the defects in the nested 
network. In the thin wall approximation, the potential 
Vj, n ((f>, x) allows the formation of three-junctions as reac- 
tions that occur exothermically, and the nested array of 
thin wall configurations is stable. In Fig. 4 we depict the 
hexagonal network of defects inside the domain wall, in 
the thin wall approximation; the dashed lines show equi- 
lateral triangles, that belong to the dual lattice. Both 
the hexagonal network and the dual triagular network 
are composed of equilateral polygons, a fact that follows 
in accordance with the Z3 symmetry. 




FIG. 4. The equilateral hexagonal network of defects, that 
may live inside the host domain wall. The dashed lines show 
the dual lattice, formed by equilateral triangles. 

We now explore the breaking of the Z2 x Z3 symmetry 
of the model. The simplest case refers to the breaking 
of the Z3 symmetry, without breaking the remaining Z2 
symmetry. We consider the case of breaking the internal 
Z3 symmetry in the following way. We take for instance 
the vacuum state vf 1 — (3/2) (1,0), and change its po- 
sition to a location farther from or closer to the other 
minima of the system, increasing or decreasing the angle 
between two of the three defects; see Fig. 5. We can do 
this with the inclusion in the potential of another term, 
proportional to the second-order power on (f>. We notice 
that the energy of the defect depends on the distance 
between the two minima the defect connects, and goes 
with the cube of it. Thus, if the vacuum state deviates 
significantly from its i^-symmetric position, we cannot 
neglect the correction to the energy of the defects. This 
changes the regular hexagonal pattern of Fig. 4 to two 
other hexagonal patterns, composed of thicker or thinner 
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hexagons. We recall that hexagonal patterns may ap- 
pear in chemical systems [ j55f , and in fluid convection fug ] 
where they may also involve non-equilateral hexagons. 




FIG. 5. The vacuum states (black dots) and the junction 
that forms the nested network, (a) and (b) illustrate the only 
two ways of breaking the Zz — » Zi symmetry. 

We explore the presence of local defects in the hexag- 
onal network by introducing penta-hepta pair of cells, 
in a local deformation of the network that disorganize 
its otherwise regular pattern. The mechanism is similar 
to that of Refs. |57||3^]. However, if the Z 2 symmetry 
that governs the host domain wall is effective, local de- 
formations may only appear in a flat surface, requiring 
the pentagons and heptagons are not regular polygons. 
This possibility may be seen in Benard-Marangoni con- 
vection; see Ref. p9| ] for a report on the experimental ob- 
servation of such paterns. But if together with the slight 
breaking of the Z% symmetry of the internal network, 
one slightly breaks the Z2 symmetry of the host domain 
wall locally, this will ultimately favor the appearence of 
local deformations composed of pair of equilateral pen- 
tagons and heptagons. Since in the network of equilateral 
hexagons, the presence of equilateral pentagons and hep- 
tagons introduce local curvature, positive and negative, 
respectively, we can understand these local defects as a 
mechanism for roughening the planar surface that con- 
tains the network. To break the symmetry of the nested 
network, one requires a slight change of position of one of 
the three minima of the nested system, so we can neglect 
the difference in energy and consider the tension as in the 
regular hexagonal network. We see that the roughening 
springs to generate higher energy states from the planar 
regular hexagonal structure. 

We now concentrate on breaking the Z 2 symmetry of 
the host domain wall. We can do this with the inclu- 
sion in the potential of a term odd in c, that slightly 
removes the degeneracy of the two minima a = ±3/2. 
Thus, the host domain wall bends trying to involve the 
local minimum, the false vacuum. To stabilize the non- 
topological structure we include charged fields into the 
system. The way one couples the charged fields is not 
unique, but if we choose to add fermions, we can couple 
them to the a field in a way such that the projection 
with {4>, x) ~ > (0, 0) may leave the model supersymmet- 
ric. This is obtained with the superpotential W(a), with 
the Yukawa coupling cPW/dcr 2 = (4/3)v / 3cr. In this case 
massless fermions bindj7Q] to the host domain wall, and 
contribute to stabilize H56[ the non-topological defect that 



emerges with the breaking of the Z2 symmetry. 

The breaking of the Z 2 symmetry can be done break- 
ing or not the remaining Z3 symmetry of the model. We 
examine these two possibilities supposing that the host 
domain wall bends under the assumption of spherical 
symmetry, becoming a non-topological defect with the 
standard spherical shape. This is the minimal surface 
of genus zero, and according to the Euler theorem we 
can only tile the spherical surface with three-junctions 
as a regular polygonal network in the three different 
ways: with 4 triangles, or 6 squares, or yet 12 pentagons. 
They are the tetrahedron, cube, and dodecahedron, re- 
spectively. They are three of the five different ways of 
tiling the sphere with regular polygons, known as the 
Platonic solids |fn}| . These three cases preserve the Z3 
symmetry of the original network, locally, at the three- 
junction points. However, if locally one slightly breaks 
the Z3 symmetry of the network to the Z2 one, the three- 
junctions can now tile the spherical surface with 12 pen- 
tagons and 20 hexagons. We think of breaking the Z3 
symmetry minimally, to the Z 2 symmetry, through the 
same mechanism presented in Fig. 5. Thus, if the symme- 
try is broken slightly we can consider the defect tensions 
as in the regular hexagonal network. 

The tiling with 12 pentagons and 20 hexagons gener- 
ates a spherical structure that resembles the fullerene, 
the buckyball composed of sixty carbon atoms 
This is the truncated icosahedron, one among thirteen 
different possibilities of tiling the sphere with regular 
polygons of two or more distinct types, known as the 
Archimedean solids ]7l{ |. The truncated icosahedron is 
one of the seven Archimedean solids constructed with 
triple junctions, and it is the one that breaks the Z3 sym- 
metry very slightly. We visualize the symmetries involved 
in the spherical structures thinking of the corresponding 
dual lattices, which are triangular lattices, but in the 
three first cases the triangles are equilateral, while in the 
fourth case they are isosceles. We recall that regular 
heptagons introduce negative curvature, so they cannot 
appear when the genus zero surface is minimal. However, 
they may for instance spring to generate higher energy 
states from the fullerene-like structure, locally roughen- 
ing the otherwise smooth spherical surface. 

We write the energy of the non-topological structure 
as E™ t = E^ t + E n , where E* t stands for the energy of 
the standard non-topological defect, and E n is the por- 
tion due to the nested network. We use E^ t = E q + Eh , 
which shows the contributions of the charged fields and of 
the host domain wall, respectively. We have Eh = Sth, 
and E n — N dt n , where S is the area of the spherical 
surface, and N and d are the number and length of seg- 
ments in the nested network. We introduce the ratio 
E^ t /E s nt = l + [N/(l+r)](t n /t h )(d/S), with r = E q /E h . 
The non-topological structure nests a network of defects, 
which modifies the scenario one gets with the standard 
domain wall. The modification depends on the way 
one couples charged bosons and fermions to the a, <$>, 
and x fields. However, if the Z% symmetry is locally 
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broken to the Zi one, the most probable defect corre- 
sponds to the fullerene or buckyball structure. But if 
the Z% symmetry is locally effective, there may be three 
equilateral structures, the most probable arising as fol- 
lows. We consider the simpler case of plane polygonal 
structures, identifying the tetrahedron (i = 3), cube 
(i = 4), and dodecahedron (i = 5). We introduce 
as the energy ratio for the i and j structures. We get 
nj = (l+r+t n /hit h )/(l+r+t n /hjt h ), fori,j = 3,4,5. 
Here h^,hi, and /15 stand for the radius of the incircle 
of the triangle, square, and pentagon, respectively. En- 
ergy favors the triangular lattice as the nested network. 
This configuration is self-dual, because the network and 
its dual are the very same triangular lattice. The two 
other configurations (that add to give the five Platonic 
solids) are the octahedron, dual to the cube, and the 
icosahedron, dual to the dodecahedron. They do not ap- 
pear in the Z2 x Z3 model because they require four- and 
five-j unctions, respectively. 

Our work can be extended in several directions. For 
instance, we could use the Z2 x Zk symmetry (k = 4, 5, 6), 
getting to /c-junctions. This allows to tile the plane with 
squares (k — 4), or triangles (k — 6), and the spherical 
surface with triangles, as the octahedron (k = 4) or the 
icosahedron (k = 5). This direction seems appropriate to 
model the recent experimental observations of squares in 
specific Rayleigh-Benard and Benard-Marangoni convec- 
tions ||74l , fr5|| . Also, in the Z 2 x Z 3 model, if the host do- 
main wall bends cylindrically, one may get to nanotube- 
like configurations [ [Z6|j77|] . As one knows, in certain types 
of nanotubes one can find polarons (78) , so we are using 
this idea to investigate the presence of chiral polarons in 
chiral nanotubes. Another line of investigations concern 
the use of real scalar fields to map laser beams inter- 
acting inside nonlinear fiber optics, as for instance in the 
model we have already presented in Ref . (79) , where laser 
beams that generate black solitons are used to give rise 
to a vector soliton of the black and bright type. 

We would like to thank C.A.G. Almeida, B. Baseia, 
FA. Brito, W. Freire, L. Losano, J.M.C. Malbouisson, 
J.R.S. Nascimento, R.F. Ribeiro, M.M. Santos, and C. 
Wotzasek for invaluable discussions, and CAPES, CNPq, 
PROCAD and PRONEX for financial support. 
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